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Outcoupling from a Bose-Einstein condensate with squeezed light to produce 

entangled atom laser beams. 

S.A. Haine and J.J. Hope 
Australian Centre for Quantum- Atom Optics, The Australian National University, Canberra, 0200, Australia^ 

We examine the properties of an atom laser produced by outcoupling from a Bose-Einstein con- 
densate with squeezed light. We model the multimode dynamics of the output field and show that 
a significant amount of squeezing can be transfered from an optical mode to a propagating atom 
laser beam. We use this to demonstrate that two-mode squeezing can be used to produce twin atom 
laser beams with continuous variable entanglement in amplitude and phase. 
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I. INTRODUCTION 

The experimental demonstration of Bose-Einstein con- 
densates (BEC) pj has lead to the development of atom 
lasers by outcoupling atoms from trapped BECs by ei- 
ther a radio frequency transition or a Raman transition to 
change the internal state of the atom to one that is either 

untrapped or anti-trapped [lHSI!HISEH[l3]. At ° m 
lasers are coherent matter waves with spectral fluxes 
many orders of magnitude higher than thermal sources of 
atoms. The coherence of these sources will enable an in- 
crease in the sensitivity of interferometric measurements 
Although current experiments usually operate in 
parameter regimes limited by technical noise, the funda- 
mental limit on these measurements will be caused by 
the shot noise of the atomic field, which will be intrin- 
sic to all interferometers without a non-classical atomic 
source. Sensitivity is increased in optical interferome- 
try by 'squeezing' the quantum state of the optical field, 
where the quantum fluctuations in one quadrature are 
reduced compared to a coherent state, while the fluctu- 
ations in the conjugate quadrature are increased. In the 
context of atom optics, it is interesting to ask whether 
highly squeezed atom optical sources can be produced. 
There is also great interest in the production of entan- 
gled atomic beams for quantum information processing 
and tests of quantum mechanics with massive particles 
[l2| . This paper will describe methods of coupling the 
quantum statistics from optical fields to produce non- 
classical atomic sources with high efficiency. 

Generation of squeezed atomic beams has been pro- 
posed by either utilising the nonlinear atomic interactions 
to create correlated pairs of atoms via either molecular 
down conversion or spin exchange collisions [l3L IT3 . fl5| , 
or by transferring the quantum state of a squeezed op- 
tical field to the atomic beam 0, 0, 0] . In the first 
case, it was shown that collisions between two conden- 
sate atoms in the \Mp = 0) state can produce one atom 
in the \Mp = +1) and one atom in the \Mp = — 1), with 
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sufficient kinetic energy to escape the trap 0, 0] . It 
was shown that this scheme produced pairs of atoms en- 
tangled in atomic spin. In the second scheme a BEC of 
molecules composed of bosonic atoms is disassociated to 
produce twin atomic beams, analogous to optical down 
conversion |15| . It was shown that the beams were entan- 
gled in the sense that phase and amplitude measurements 
on one beam could infer phase and amplitude measure- 
ments of the other beam better than the Heisenberg limit. 
Although each atomic pair is perfectly correlated in direc- 
tion in each of these schemes, there is very little control 
of direction of each pair, so the spectral flux would be 
limited. 

The generation of nonclassical light is well established 
experimentally [l6|. This suggests that a nonclassical 
atom laser output could be generated by transferring a 
the quantum state of an optical mode to an atomic beam. 
Moore et al. showed that a quantized probe field could 
be partially transferred to momentum 'side modes' of a 
condensate consisting of three-level atoms in the presence 
of a strong pump field |17| . Jing et al. performed a single 
mode analysis of the atom laser outcoupling process for 
a two-level atom interacting with a quantized light field, 
and showed that the squeezing in light field would oscil- 
late between the light field and the atomic field at the 
Rabi frequency 0|- As this was a single mode analysis, 
the interaction with the atoms as they left the outcou- 
pling regi on was not taken into account. Fleischhauer 
et al. 19] showed that Raman adiabatic transfer can be 
used to transfer the quantum statistics of a propagating 
light field to a continuously propagating beam of atoms 
by creating a polariton with a spatially dependent mix- 
ing angle, such that the output contained the state of the 
probe beam. 

In this paper, we model the dynamics of an atom laser 
produced by outcoupling three-level atoms from a BEC 
via a Raman transition, and investigate the transfer of 
quantum statistics from one of the optical modes to the 
atomic field. Ideal transfer will occur when the time 
taken for each atom to leave the outcoupling region is 
a quarter of a Rabi period. The finite momentum spread 
of a trapped condensate means that there will be a broad- 
ening of the time taken to leave the outcoupling region, 
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and hence ideal transfer will not be possible. To deter- 
mine the effectiveness of the quantum state transfer, we 
require a multimode model that takes into account back 
coupling and the finite momentum spread of the conden- 
sate. 

In Section II we describe an atom laser beam made by 
outcoupling from a BEC using a non-trivial optical mode, 
using the simplest possible model that contains the spa- 
tial effects in the output mode. We derive the Heisen- 
berg equations of motion for this system under suitable 
approximations. Section III introduces the method used 
to solve these equations and investigates some properties 
of the outcoupled atoms, showing that complicated spa- 
tial behaviour occurs in the output even when the optical 
and BEC fields are described by a single mode. In sec- 
tion IV we investigate continuous outcoupling with two 
mode squeezing, and show that it can be used to gener- 
ate continuous variable entanglement in twin atom laser 
beams propagating in different directions. 




(trapped) (untrapped) 

FIG. 1: Internal energy levels of a three level atom. A con- 
densate of state |1) atoms confined in a trapping potential 
are coupled to free space via a Raman transition affected by 
a probe beam (annihilation operator &13) which is detuned 
from the excited state (|3)) by an amount Ai, and a pump 
field (annihilation operator 023) which is detuned from the 
excited state by an amount A2. 
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II. OUTCOUPLING USING A NONCLASSICAL 
OPTICAL FIELD 



When an atomic and an optical held are coupled, and 
they can both be described by a single mode, then com- 
plete state transfer must occur between them in a Rabi- 
like cycle. When producing an atom laser beam in this 
manner, however, the single mode approximation cannot 
be made for the output held, even though it may be ap- 
plicable to the optical and BEC fields. In this section we 
develop such a model, and derive the Heisenberg equa- 
tions of motion for the output held and the optical held 
operators. 

We model an atom laser in one dimension as a BEC 
of three-level atoms coupled to free space via a Raman 
transition, as shown in figure^ State |1) represents the 
internal state of the trapped condensate, |3) the excited 
state, and |2) the untrapped atomic mode. 013 is the 
annihilation operator for the probe optical mode (tran- 
sition |1) — > |3)), and 023 is the annihilation operator 
for the pump optical mode (transition |2) — * |3)). The 
pump held is assumed to be a large coherent state, much 
stronger than the probe held, so it is approximated well 
by a classical held 323023 = ^23 e ~ i( "~ A2)t - Tne Hamilto- 
nian for the system (in the rotating wave approximation) 



7~L 7~tatom flight ~t~ 'H-atom— light 

2m 

h 2 k 2 

+ / + tuj)i>3(k)dk + h{uj - Ai)a f 13 ai3 



(1) 



$(fc)#o^i(fc)* + ^- / k 2 4>l(k)4> 2 {k)dk 



+ ft.9i3 J ipi(k)ipl(k + k 13 )d 13 + ip[(k)ip 3 (k + k 13 )a[ 3 dk 

+ h ( n 23 ^ 2 (k)Mk + k 23 )e l(ul - A2)t 



where tpi(k) is the k-space annihilation operator the con- 
densate mode (internal state |I)), ^ 3 {k) is the annihila- 
tion operator for atoms in the excited atomic state (|3)), 
and ip2{k) is the annihilation operator for the untrapped 
free propagating mode (|2)). The annihilation operators 
obey the usual bosonic commutation relations: 

[&(*), = $(#)] = 0, (2) 

Hq is the single particle Hamiltonian for the trapped 
atoms, m is the mass of the atoms, tt 23 is the Rabi 
frequency for the pump transition, gi 3 is the coupling 
strength between the atom and the probe held, huj is the 
internal energy of the excited state |3) atoms, and hki 3 
and hk 23 are the momentum kicks due to the pump and 
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probe light fields respectively. For simplicity we have as- 
sumed a laser geometry where the pump and probe fields 
are counter propagating, to give the maximum possible 
momentum kick to the untrapped atoms. The equations 
of motion for the Heisenberg operators are: 





= -r- ^i(fc) + gi>3(k + kiz)a) 
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r0 2 (fc) 
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= -Z—i>2(k) + Q23-4>3(k + k 23 ) 

Am 


(4) 


iijj 3 (k) 


hk 2 

= {— + A 2 )^ 3 (k) + g 13 Mk-ki 3 )a 
Am 

+ ^ 23 ip 2 (k ~ k 23 ) 
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ia 


= 6a + 313 J ' i>l[k - ki 3 )Tp 3 {k)dk 


(6) 



Where ijj 3 = t/5 3 e*( u ' -Aa J* and a = a^e^-^* , and <5 = 
(A 2 — Ai) is the two-photon detuning. 

The population of state |3) will be much less than 
the other levels when the detunings (Ai, A2) are much 
larger than the other terms in the system (including 
the kinetic energy of the excited state atoms). Further- 
more, most of the dynamics will occur on time-scales 
less than so in this regime we can set tp 3 (k,t) w 

^7(513^1 (fc -k 13 ,t)d + f^ 3 V) 2 (fc -ki3,t)). If the conden- 
sate has a large number of atoms and is approximately in 
a coherent state, we can write ij)i{k,t) « y/N(f>o(k)e~ Wtt , 
where <po(k) is the condensate wavefunction (which we 
will assume is in the ground state of the harmonic oscil- 
lator, with uj t the trapping frequency) and N is the con- 
densate number. We have ignored the atom-atom inter- 
actions in our model, which is valid only if the condensate 
is dilute. Strong atom-atom interactions would have the 
effect of introducing complicated evolution to the quan- 
tum state of the condensate mode. Inclusion of these 
effects is not possible with our method, and a more com- 
plicated technique such as a phase space method would 
be required 20] . The approximation of ignoring the back 
action on the condensate is only valid if we are in the 
regime where the outcoupling is weak, ie. the number 
of photons in the probe field is is much less than the 
number of atoms in the condensate. In an experiment, 
measuring the quantum noise on the atom laser beam 
would require small classical noise on the beam, and in 
practice this is easier to achieve with weak outcoupling. 
With these approximations our equations of motion for 
the free propagating atoms and the probe field become 

#(fc) = Lu {k)tlj(k) - Q (k)a (7) 
ia = uj a a — J ilg(fc)-0(fc)c?fc (8) 

with 4,{k) = Mk)e^ 1 , Mk) = - ^ - wt), 

"a = )' and Q o{k) ^ gi 3 VN^Mk+k 23 -k 13 ). 

In the next section we will discuss the solution to these 
equations and the properties of the outcoupled atoms. 



III. PROPERTIES OF THE OUTCOUPLED 
ATOMS 

The solution to equations J7J and JSJ) is 

i>(k,t) = J f(k,k',t)Mk)dk' + g(k,t)a a (9) 
a(t) = p(t)a s + J q{k\t)i> a {k')dk' (10) 

Where a s = a(t = 0) and ip s (k) = i/j(k,t = 0) are 
the Schrodinger picture operators, and f(k, k', t), g(k, t), 
p(t), q(k',t) are complex functions satisfying: 

if(k,k') = uj Q (k)f(k,k')-n (k)q(k') (11) 
ig(k) = u (k)g(k) - n (k)p 

ip = u> a p — J fl'Q(k)g(k)dk 
iq(k') = u a q{k')- [ {l* (k)f(k,k')dk 



with initial conditions f(k, k',t = 0) = 6(k — k'), p(t = 
0) = 1, and g{k,t = 0) = q(k',t = 0) = 0. This ansatz 
has reduced the field operator equations to a set of cou- 
pled partial differential equations. This will only be pos- 
sible for Heisenberg equations of motion that do not have 
terms with products of operators, but it allows the pos- 
sibility of an analytic or numerical solutions to the full 
quantum problem. 

We solved equations (|llfl numerically using a fourth 
order Rungc Kutta algorithm using the XMDS numeri- 
cal pacakge |2J . We chose parameters realistic to atoms 
optics experiments with Rb 87 atoms. Unless stated oth- 
erwise, we have set m = 1.4x 10~ 25 kg, cu t = 0.25 rad s _1 , 
|k23 — ki3 1 = 1.6 x 10 7 m _1 , which corresponds to twice 

the wave number of the 2 Sj. ,F= 2 — * 2 P 3 ,F= 3 transition 

2 2 

in Rb 87 . <j>o(k) was chosen to be the (normalized) ground 
state momentum space wave function of a condensate 
(ignoring interactions) in a harmonic trap, and we set 
fio(fe) = tt(t) Q (k-k2 3 -ki 3 ) with ft = 90rads _1 . There- 
suits are reasonably insensitive to the absolute magnitude 
of Lu a and ujq, but they are quite sensitive to the relative 
values. To maintain resonance between the two fields, we 

set ^25J~ = h ^ k23 2m 1 ^ — ^a—^t- These relationships can 
be obtained with physically realistic parameters and are 
consistent with the approximations made in this model. 
We set uj a = 20 rad s _1 . 

Figures © , 10 , and (J5J show the solutions to equa- 
tions (fTTf) for the values indicated above. 

The solution of equations (| 1 1 1> gives us the solution of 
Eqs.© and (|10fl for all possible initial quantum states of 
the optical field and the free propagating atomic field. 
We will assume that the initial state of the field is 
1-0) = \light) ® {|0)}fe , where \light) represents an ar- 
bitrary state for the optical mode, and {|0)}/- represents 
a vacuum mode at all points in k space for the atomic 
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FIG. 2: \f(k,k',t = 0. lis) | 2 for the values of parameters 
indicated in the text, kkick is the kick acquired due to the 
Raman transition, ie kkick = &23 — fci3- The function was 
discretized for numerical calculation by replacing S(k — k') 

with ^ , where Afc is the grid spacing. The dip in the 
function near the k = k^ick resonance shows how the quantum 
state of the atoms has been affected by the interaction with 
the optical fields. 



FIG. 4: |p(t)| 2 as found numerically for the values of param- 
eters indicated in the text. 





FIG. 3: \g(k,t)\ 2 as found numerically for the values of pa- 
rameters indicated in the text. This shows that atoms are 
created around k = kkick with a quantum state related to the 
initial state of the probe field. 



FIG. 5: \q(k',t)\ 2 as found numerically for the values of pa- 
rameters indicated in the text. 



field. The expectation value of the density of outcoupled 
atoms (* t (a;)* (x)) with = ^= / iji(k)e ikx dk is 

(¥(x)*(x)) = \G(x)\ 2 (ata s ), G(x) = ^J g{k)e lkx dk 

(12) 

It is interesting to note that when the initial state of the 
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FIG. 6: (* t (a;)*(s)) for = 1000. 



untrapped atomic field is the vacuum, then the spatial 
structure of the density of the untrapped atoms at later 
times depends only on the the functional form of G(x), 
which depends on the efficiency of the outcoupling pro- 
cess. Figure © show the density of outcoupled atoms 
when the expectation value of the initial number of pho- 
tons is (al& 8 ) = 1000. 

The number operator is N = J (x)^(x)dx. Using 
our solution (Eq. and our initial state, the expectation 
value is (N) = (a\a s ) J \G{x)\ 2 dx. The variance of the 
number operator is 



V(N) 



(N 2 ) - (A>) 2 

(¥(x')$(x')&(x)$(x))dxdx' 

(&(x')¥(x)$(x')i>)(x)dxdx' 



)dx 



(13) 

Using our solution (Eq. EJ) and our initial state, this 
becomes 



V(N) = N*((ala%a s ) - (a%) 2 ) + Nafta.) 

= N 2 V (aids) + N G (1-N G ) (ata s ), (14) 



with N G = J \G(x)\ 2 dx. We note that as N G -> 1, the 
variance in the number of outcoupled atoms approaches 
the variance of the initial optical mode, as the quantum 



FIG. 7: The relative variance v(N) 



for the out- 



(JV> 

coupled atoms versus time for different states of the optical 
field. The solid line represents the initial optical mode in a 
coherent state \a) with \a\ 2 — 1000, the dotted line represents 



a squeezed state \a,r) with |a| 



1000, 



1.38, and the 



dashed line represents a Fock state \n) with n = 1000. 



statistics of the outcoupled atoms depends only on the 
initial quantum state of the optical held and the efficiency 
of the outcoupling process. Figure (7J shows the variance 
of the outcoupled atoms versus time for different states 
of the optical mode. A more interesting observable to 
look at is the flux of the outcoupled atoms, as the spatial 
structure of the outcoupled beam becomes apparent. The 
flux operator is: 



J(x) 



ih 

2m 



V* t (x)#(x) - &(x)V&(x] 



(15) 



which, using our solution for i/)(k) becomes 



J(x) 



J f (x, k', k")'4>l(k')i> s (k'')dk'dk'' 
+ J g (x)dld s 

Jf g (x,k')Tpt{k')d s dk' 

J gf (x,k")i> s (k")aidk" 



(16) 



with 
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The variance in the flux has two terms, one propor- 
tional to the variance in the photon number, and the 
other proportional to the photon number itself. For 
a Fock state photonic input the first of those terms is 
zero, and the variance is proportional to the function 
J Jgf(x,k)Jf g (x,k)dk. This can be contrasted to the 
case where the optical field is in a coherent state, and 
the total variance in the flux is simply proportional to 
the function J 2 + J J g f(x,k)Jf g (x,k)dk. A reasonable 
measure of the transfer of the quantum state of the zero- 
dimensional photon field to the larger space of the output 
pulse is therefore the function 



v(J) 



I Jgf(x,k')J fg (x,k')dk' 
J g+ I J gf( x ' k ) J fgi x i k ) dk ' 



(20) 
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FIG. 8: Flux of outcoupled atoms at a point in the atomic 
beam (x — 1.5 mm) for (aja s ) = 1000. 



J f (x, k', k") = ^- (yF*{x, k')F(x, k") - F* (x, k')VF(x, k" 
J„(x) = ^—(wG*(x)G(x) - G*(xWG(x)) 
J fg (x, k') = ^-(vF*(x, k')G{x) - F*(x, fe')VG(^)) 
J gf (x, k") = ^ (VG* (x)F(x, k') - G*(x)VF(x, k")) 

and with F(x, k') = J f(k, k')e tkx dk. 

Using our initial state \ip) = \light) (g) {|0)}fc, the ex- 
pectation value of the flux operator becomes: 



(J(x)} = Jg(x)(ala 



(17) 



This shows that the mean atom flux in the output pulse 
depends only on the details of the coupling process, and 
not on the statistics of the outcoupling field. Figure 03 
shows the flux of outcoupled atoms for (aja s ) = 1000. 

To investigate how the quantum statistics are trans- 
ferred to the atomic beam, we look at the variance in the 
flux. 



V(J) = (J 2 ) {Jf 



Jg(ala s ala s ) - J 2 g {d\d 8 ) 2 



(18) 

+ J J Jg f (x,k')J f g(x,k")(atMklaJ4(k''))dk'dk'' 
J 2 g V{ala s ) + (ala s ) J J gf (x,k')J fg (x,k')dk' (19) 



which shows the minimum possible variance in the output 
flux normalised to the flux variance produced by output 
with a coherent optical state. 

Figure (JHJl shows v (J) for different values of the cou- 
pling constant Q. Even in our simplified model where 
we have assumed a single mode for the optical beam 
and the condensate, the outcoupled atoms still display 
complicated spatio-temporal dynamics. Weak outcou- 
pling gives a steady flux, but very little suppression of 
the shot noise because the timing of the output of each 
atom becomes uncertain, making the number statistics 
uncertain in the transient period. When the outcoupling 
j^te is increased, a significant amount of flux squeezing 
is displayed in a localised pulse. Further increase of the 
outcoupling rate shows more complicated dynamics, as 
some of the outcoupled atoms are coupled back into the 
condensate. This causes the atoms to come out in a 
series of pulses, with less flux squeezing than for opti- 
mal outcoupling. An interesting sidenote is that the flux 
variance produced by the coherent optical state (the de- 
nominator of v(J)) is simply proportional to the flux it- 
self, with the same proportionality constant for all times, 
and all values of Q. Although the variance in the num- 
ber of the pulse is quite insensitive to the strength of 
the coupling between the trapped and untrapped fields, 
this is not true for the variance in the flux. Figure l|10|) 
shows the maximum suppression of shot noise in v(J) 
for different values of f2. We can estimate the outcou- 
pling that will produce the minimum v( J) by finding the 
maximum Rabi frequency that will not cause significant 
back-coupling to the condensate. Equating the quarter- 
period of a Rabi oscillation Tuabi/^ = 7r/ (2S1) , with the 
time taken for the kicked atoms to leave the coupling 

region T leave = ^f^ m /(h\k 23 - k 13 |) where 

is the spatial width of the condensate. From this we 
can estimate that optimum outcoupling will occur when 
fl w ^l^-k^l « 250 rad s _1 for the parameters used in 

this paper. This agrees well with the calculated minimum 
shown in figure ( HUP . 

We have shown that the quantum statistics of an op- 
tical mode can be transferred to an atom laser beam to 
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FIG. 9: v(J) at a point in the path of the atomic beam 
(x = 1.5 mm) for SI = 18 rad s^ 1 (dotted line), Q. = 144 rad 
s" 1 (solid line), and SI = 270 rad s _1 (dashed line). Coupling 
weakly produces a long pulse, but the variance in the flux is 
almost unaffected by the statistics of the optical state. Cou- 
pling too strongly causes significant back coupling from the 
output field to the photonic state. 




FIG. 11: Twin atom laser beams produced by outcoupling 
with two- mode squeezed light. An OPO is driven by a clas- 
sical, non-depletable driving field (/3). The x' 2 ' process pro- 
duces two optical modes ai and which are used to outcou- 
ple the atom laser beams. 
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FIG. 10: Minimum value of v( J) versus Q . The shot noise is 
below the vacuum noise for all values of £1 when using a Fock 
state to outcouple. 



produce a pulse of atoms with better defined number, or 
to partially suppress the fluctuations in the flux. Fur- 
thermore, we have shown that the quantum statistics of 
the optical mode can be transferred independent of the 
initial quantum state of the optical mode, which suggests 
that two-mode optical squeezing could be used to gener- 
ate spatially separated entangled atomic beams. In the 
next section we investigate the possibility of using twin 
optical beams produced from a non-degenerate optical 
parametric oscillator to generate two entangled atomic 
beams propagating in different directions. 



IV. EPR BEAMS 

Continuous wave generation of correlated atom beams 
requires a more complicated scheme. We consider a probe 
field created by a nondegenerate OPO producing twin 
optical beams (figure l|llfl ). These modes have the same 
wavelength, but travel in different directions and hence 
they have different momenta. The OPO is driven by a 
classical, non-depletable driving field. This will produce 
twin atom laser beams with different momenta. This 
differs from the previous case in that it allows continuous 
outcoupling of the atoms, rather than just a pulse. The 
Hamiltonian for the system is now 
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W = TLatom (21) 

+ H(oj — Ai)a\ai + h(oj — Ai)a 2 &2 

+ +/3*aia 2 e^*) 

+ %13 / ^i(ft)V>l(ft + fti)ai + ip\(k)il> 3 (k + kt)a\dk 



+ %13 / i>i(k)ipt{k + k 2 )d 2 + i>l(k)ip 3 (k + k 2 )aldk 

+ h J n 23 ^ 2 (k)Mk + fcoK ( "- A2)t 

+ r223V>2(fc)^(fc + fco)e _i( ^ A2) *dfc 



where ai and S2 represent the annihilation operators for 
the twin probe fields produced from the down conversion 
process, both assumed to affect the |1) — > |3) transition, 
(3 is the complex amplitude of the pump field, lo p is the 
frequency of the pump, and x is the nonlinear coeffi- 
cient of the down conversion medium. Kk\ and Hk 2 are 
the magnitudes of the momentum kicks due to absorp- 
tion from photons in d\ and d 2 respectively. We have 
assumed that the photons are resonant in an optical res- 
onator with 100% reflective mirrors. This assumption is 
valid as the dominant form of loss out of the cavity will 
be due to atomic absorption. For computational con- 
venience in our one-dimensional model, we have chosen 
ki — ko = — (k2 — ko), ie. the resultant momentum kicks 
that the atoms obtain after being outcouplcd are of equal 
magnitude and opposite direction. By adiabatically elim- 
inating the excited state, and assuming the condensate is 
a large coherent state as before, we obtain the following 
equations of motion for the outcoupled atoms and the 
probe fields: 



ii/)(k) 
iai 

id 2 



with ip(k) 



w (k)i>(k) - Qi(fc)ai + Cl 2 (k)a 2 (22) 
w a ax - f Sl*(k)i>(k)dk (23) 

u a a 2 - / £l* 2 (k)i;(k)dk (24) 



%(fc) 



+ xl3a\ et(2{uJ ^ A2) ^ p)t - n * c ai 



i(u-A a )t 



S 2 A 



^ 2 (k)e %{ , Uj 

\- ko — kj) for j = 1,2, and flc = 

~aT I ^oC — k\)4>(j{k — k 2 )dk. The Qc cross coupling 
term between the two optical modes is due to atoms ab- 
sorbing a photon from one beam and emitting it into the 
other beam. This term will be small due to the large 
momentum difference between the two modes. However, 
the functional form of Qc is due to our assumption that 
the condensate remains single mode. Cross coupling be- 
tween the two optical modes will cause momentum 'side 
bands' on the condensate mode ,17], but the effect of 
this cross coupling will be small if the number of photons 
in the probe beam is small compared to the number of 



atoms in the condensate. As the results in this section 
are calculated in a parameter regime where the chance of 
an outcoupled atom coupling back into the condensate is 
low, it is valid to neglect this term in our calculations. 
The general solution to equations (1221) is 



(25) 



= / f+{k,k\t)i) s (k)dk' 

/_ (ft, ft', t)$ (ft)dfc' + gi+ (k, t)a ls 

+ 51- (ft, t)a\ s + g 2+ (ft, t)a 2s + g 2 - (ft, t)a\ s 
di(t) = pi+(t)au +pi-(t)a\ s 
+ P2+{t)a 2s +p 2 -(f)4« 

+ / P3+ (k' \t)i> a (k')dk' - 

a 2 (t) = q 1+ (t)a ls + q 1 -(t)d\ s 
+ q 2 +(t)a 2s +q 2 -(t)al s 

q 3+ (k',t)Mk')dk'- 



(26) 



p 3 -(k' ,t)^l(k')dk' 
(27) 



g 3 -(fc' ,t)i>t(k')dk' 



where f±(k,k',t), 9i, 2 ±(k,t), Pi, 2 ±(t), p 3 ±(k',t), 
Qi,2±{t), q 3 ±(k / ,t) are complex functions satisfying dif- 
ferential equations obtained by substituting the solutions 
(|25|l into (|22|l (see appendix) . From the solution of these 
equations, we can calculate any observable of the system. 

We solved the equations l|34[l numerically for xft — 
80s" 1 , rij(k) = Q,(f)o(k -k - kj) with ft = 108 rad s" 1 
for j = 1,2. We set ki~k = -(k 2 -k ) = 1.6x 10 7 m -1 , 
and ujp — 2(cj — A2) = 2uj a , with all other parameters as 
before. If we assume that the two optical modes and the 
untrapped atomic field are initially in the vacuum state, 
using equation the expectation value of the atomic 
density p[x) = (* 1 \x)) is 

p{x) = (&(x)^(x)) (28) 
\F_(x,k')\ 2 dk' + \G^(x)\ 2 + \G 2 -(x)\ 2 



Where F_{x,k') = -j= J f^{k,k')e lkx dk, G x _(x) 



-2=f gi „(ky k *dk and G 2 .{x) = J= f g 2 _( k )e ik *dk. 



Figure l|12|l shows the atomic density versus time. Two 
atomic beams in opposite directions are produced, with 
steady flux. To check whether there are correlations 
present in the two atomic beams, the relevant observ- 
able is the difference in the flux of the two beams. If the 
two beams are completely uncorrelated, then we would 
expect 



V(J(x ) - J(-xo)) > 2V(J(x )) 



(29) 



Where xq and —xo are points that lie in the rightward 
and leftward propagating beams respectively. Figure 
(|13J) shows the variance of the flux difference V(J(xo) — 
J(—xo)) versus time, and shows that the fluctuations in 
the difference are approximately eight times less than for 
uncorrelated atoms. 
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the quantum limit. We define four quadratures: 

X± = J(L* ± (x)^(x)+L ± (x)¥(x))dx (30) 

Y± = i J(L* ± (x)$(x) - L±{x)&(x))dx (31) 

(32) 
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FIG. 12: Density of outcoupled atoms versus time. Outcou- 
pling using light from a non-degenerate OPO produces two 
atomic beams in opposite directions, with steady flux. 
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FIG. 13: Variance in the flux difference. The fluctuations are 
eight times smaller than for uncorrelated atoms. 



We now investigate whether we can use this process to 
generate entanglement between the two atomic beams. 
We quantify our entanglement using the EPR criterion 
of Reid and Drummond p2l | , the requirement being that 
two conjugate variables on one of the beams can be in- 
ferred from measurements on the other beam to below 



L±(x,t) 



y/\xi - X 2 \ ' 

= otherwise 



if ± X\ > X > ±X2 



(33) 



The commutator of the conjugate quadratures give 
us the uncertainty relation V(X±)V(Y±) > 1 since 
/ \L±(x)\ 2 dx — 1. The beams are entangled under the 
EPR criterion if by making measurements of quadra- 
tures of one beam (eg. X+ and then quadra- 
tures of the other beam (Jf_ and YL) can be in- 
ferred to better than this quantum limit. Quantita- 
tively: V in f(X-)V in f(Y_) < 1 is the requirement for 

entanglement, where V mf (X ± ) = V(X ± ) - (y( **/f ))2 , 



V inf (Y±) = V{Y±) 



and V(a, b) = {ab} 



V(X T ) 

(a) (b) . We note here that the correlations present are 
between conjugate quadratures of each beam. Fig- 
ure (|14(l shows the product of the inferred variances 
V inf (X_)V inf (Y-) plotted against time. As the inten- 
sity of the beams increase and become more monochro- 
matic, the product of the inferred variances dip well be- 
low the requirement for entanglement. The initial in- 
crease is due to the beams initially not approximating 
plane waves. This could be fixed by appropriate choice 
of L±(x), to better match the mode shape of the output 
atom laser beams. In the long time limit, the product 
of the inferred variances is on the order of three orders 
of magnitude below the classical limit, demonstrating 
that this system produces an almost pure EPR correlated 
state. Our model uses an ideal OPO, so the squeezing in 
the optical modes would grow without bound in the ab- 
sence of the damping due to the atoms. In practice, the 
entanglement on the atomic beams will not exceed the 
optical entanglement that can be obtained from a real 
OPO. The limit to the entanglement between the atomic 
beams in this model is given by the finite momentum 
width of the condensate. As the EPR state is expected 
to be very pure, the dominant noise in an experiment may 
actually be due to some of the effects we have ignored in 
our model due to their small effect on the dynamics. In 
particular, there may be a small reduction in fidelity due 
to effects of the back action of the outcoupling on the 
condensate wavefunction, which we have ignored in this 
calculation. 
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P3-(k',t), qi+(t), qi-(t), q 2+ {t), q 2 -{t), q 3+ (k',t), 
q 3 -(k',t), must satisfy: 




FIG. 14: Product of the inferred variances 

V inf (X-)V inf (Y-) versus time. As the system goes to 
steady state, the requirement for entanglement is satisfied. 



V. CONCLUSION 

We have modelled the dynamics of an atom laser pro- 
duced by outcoupling from a Bose-Einstein condensate 
with squeezed light. We modelled the multimode dy- 
namics of the output field and showed that a significant 
amount of squeezing can be transferred from an optical 
mode to a propagating atom laser beam. We also demon- 
strated that two-mode squeezing can be used to produce 
twin atom laser beams with continuous variable entan- 
glement in amplitude and phase. 

This research was supported by the Australian Re- 
search Council Centre of Excellence for Quantum Atom 
Optics. We would like to acknowledge useful discussions 
with M. K. Olsen, A. M. Lance and H. A. Bachor. 



VI. APPENDIX 

f + (k,k',t), f-(k,k',t), g 1+ (k,t), gi-(k,t), g 2+ (k,t), 
g 2 -{k,t), pi+(t), pi-(t), p 2 +(t), P2-(t), p 3+ (k',t), 



if+(k,k') 
if-(k,k r ) 

m+(k) 

igi-(k) 
ig 2 -(k) 

ipi- 
iP2+ 
ip2~ 
iP3+(k') 
ip 3 _(k') 
iqi+ 
iqi- 
iq 2 + 
iq 2 - 
%+(&') 
iq 3 _(k') 



^ (fc)/+(fc,fc')-^i(fc)p 3+ (A:') (34) 
il 2 (k)q 3+ (k') 

Lu (k)f-(k, k') - n 1 (A)p 3 _(fc / ) - n 2 (k)q 3 _(k') 

u>o(k)gi+(k) - fli(k)p 1+ - il 2 {k)q 1+ 
uo(k)gi-(k) - Cli(k)p!- - il 2 (k)qi- 
uo(k)g 2 +(k) - Qi(k)p 2 + - il 2 (k)q 2+ 
u> {k)g 2 -(k) - il\{k)p 2 - - il 2 (k)q 2 - 

iv a p 1+ - J il\{k)g l+ {k)dk + XpQi- 
n*i(k)gi-(k)dk + x P 1i+ 



= LOaPl 



= U>, 



= V, 



uj a p 2 + - J il{{k)g 2+ {k)dk + x P q 2 - 
uj a p 2 - - J HI {k)g 2 - (k)dk + x P q 2 + 
Lo a p 3+ {k') - J ill(k)f+(k, k')dk + x P q* 3 -(k') 
u a P3-(k') - J tt*i(k)f-(k, k')dk + x P <ll + (k') 

a qi+~ J ^l(k)gi+(k)dk + x P Pl- 
u a qi- - J ill (k)gi- {k)dk + x P Pi+ 
uj a q 2+ - j ill(k)g 2+ (k)dk + x P P 2 - 
uj a q 2 - - J ill {k)g 2 - {k)dk + x P P 2 + 
u a q 3 +{k') - j ill(k)f+(k,k')dk + X P P 3 -(k') 
a q 3 -{k')- J ill(k)f_(k,k')dk + Xp P* 3+ (k') 



with x P = xfl^ 2 ^ ^ " p ^*! an d initial conditions 
f+(k,k',t = 0) = 5(k - k'), p 1+ (t = 0) = 1, q 2+ (t = 
0) = 1 with all other fields zero. From the solution of 
these equations, we can calculate any observable of the 
system. 
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